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X-ray and strong lensing mass estimate of MS2 137.3-2353 
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ABSTRACT 

We present new mass estimates of the cluster of galaxies MS2 137.3-2353, inferred from 
X-ray and strong lensing analyses. This cluster exhibits an outstanding strong lensing 
configuration and indicates a well-relaxed dynamical state, being most suitable for a mass 
reconstruction which combines both techniques. Despite this, several previous studies have 
claimed a significant discrepancy between the X-ray and the strong lensing mass estimates. 
The primary aim of this paper is to address and explain this mismatch. 
For this purpose, we have analysed Chandra observations to recover the profiles of the 
intra-cluster medium properties and, assuming a functional form for the matter density, 
the total mass distribution. The notable strong-lensing features of MS2137.3 allow us 
to reconstruct its projected mass in the central regions with good accuracy, by taking 
advantage of the lensing inversion code Lenstool. We compare the results obtained for both 
methods. Our mass estimates for MS2 137.3 are in agreement within errors, leading to a 
mean, extrapolated value of M200 — 4.4 ± 0.3 x 10 14 M Q , under the assumption of the 
Navarro-Frenk- White (NFW) mass profile. However, the strong lensing mass estimate is 
affected by the details of the BCG mass modeling, since the radial arc is a very sensitive 
probe of the total mass derivative in the central region. In particular, we do not find ev- 
idence for a high concentration for the NFW density profile, as reported in some earlier works. 

Key words: galaxies: clusters: individual, MS2137. 3-2353 - cosmology: dark matter, gravi- 
tational lensing - X-ray: galaxies: clusters. 



1 INTRODUCTION 



According to the hierarchical growth of cosmic structures im- 
plied by the concordance model, clusters of galaxies are the 
most recent gravitationally-rel axed s t ructures. They are ideal cos- 
mological tools (see e.g. ivoitl 120051 ; iBorganil 120061 : iTozzil |2007| ; 
Diafe rioT Schindler & DolagT 2008h . provided we have a reliable 
measument of their masses. The most promising techniques for 
estimating cluster masses are via X-ray and gravitational lensing 
analyses. 

Cluster X-ray luminosity is mainly due to diffuse plasma emis- 
sion: the deep potential wells produced by the dark matter com- 
ponent heat the intra-cluster medium (ICM) to temperatures of 
some keV, at which point X-ray photons are emitted through ther- 
mal bremsstrahlung. The X-ray emissivity is proportional to the 
square of the gas density, thus it is an excellent tracer of the three- 
dimensional cluster potential. However, X-ray measurements of 
cluster masses imply the assumption of hydrostatic equilibrium of 
the ICM with the dark matter potential and of spherical symmetry 
of the cluster mass distribution. Hence the total mass profile can be 
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inferred from the radial profiles of temperature and gas density. 
On the contrary, the gravitational lensing effect, i.e. the light deflec- 
tion due to the space-time curvature induced by massive objects, 
allows for the determination of the projected surface mass density 
of the lens, regardless of its dynamical state or the nature of the in- 
tervening matter. However, this effect is determined by all massive 
structures along the line of sight, so lensing mass measurements 
are subject to foreground and background contaminations. More- 
over, the lensing effect is very sensitive to features of the mass dis- 
tribution such as its ellipti city and asymmetry, a s well as to the 
presence of substructures (Meneghetti et al. 2007), which compli- 
cate the mass reconstruction. 

The X-ray and lensing methods are clearly complementary, al- 
lowing in principle to combine 2- and 3-dimensiona l constraints 
to obtain a reliable mass estimate (see for ex ample Allerj [l998l : 
lEttori & Lombardil 120031 : iBradac et alj l2008h : a limit to such 
a comparative analysis is the significant disagreement between 
strong l ensing and X-ray mass estimates cl ai med in the lit- 
erature JWu& Fand 1 19961: ISmail et al.l 1 1991 lOtaetalJ |2004 
IVoigt & Fabiarjl2006l : lGTtti et alJl2007l ; lHalkola et al]|2008h . Many 
convincing explanations have been suggested, but to date, the dis- 
crepancy between X-ray and strong lensing results appears to be an 
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Moreover, galaxy clusters enable one to determine dark matter 
(DM) halo density profiles, being mainly composed by dark mat- 
ter. As a consequence, they could help in verifying the CDM model 
predictions for the formation of cosmic structures: for example, if 
DM density profiles can be approximated on all mass scales by 
"uni versal" profiles, such as the Navarro- Frenk- White (NFW) pro- 
file ^Navarro. Frenk & White! 199511 19971) , as predicted by N-body 
simulations of structure formation in hierarchical collapse mod- 
els. In particular, the inner DM density slope is a very debated 
issue: several stu dies have argued that it could deviate from the 
NFW predictions dMoore et al.lll998l : |jing & Sutoll2Q00l:lGao et al.1 
2008). Galaxy clusters play a crucial role in solving this issue, 
since they allow to distinguish between cuspy and flat profiles 
due to their low concentrations. An emblematic case in both of 
these debates is MS2137. 3-2353: even though this cluster has a 
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In this work we will focus on the latter topic. We present new X- 
ray and strong lensing mass estimates, based respectively on a non- 
parametric analysis of gas density and temperature profiles through 
Chandra observations, and on a strong lensing reconstruction, per- 
formed using th e Lenstool analysis software dKneibetalJll996l ; 
ljullo et alj2007h . 

This paper is organised as follows. In §[2] we will discuss our X- 
ray analysis, focusing on data reduction and on the method applied 
to recover the total and gas mass profiles, and we will summarize 
our X-ray results. The strong lensing analysis is presented in § [3] 
where we will briefly discuss our main findings. We will compare 
the X-ray and strong lensing results in § [4] a comparison with pre- 
vious analyses can be found in §|5] Finally, we will summarize our 
results and draw our conclusions in §|6] 

Throughout this work we assume a flat ACDM cosmol- 
ogy, with matter density parameter f2 m =0.3, cosmological con- 
stant density parameter r?A=0.7, and Hubble constant Hq = 
70 h _1 km s _1 Mpc _1 . Unless otherwise stated, all uncertainties 
are referred to a 68% confidence level. 



2 X-RAY ANALYSIS 
2.1 X-ray data reduction 

MS2137.3-2353 is a rich, X-ray luminous cluster at z = 0.313 
dStocke et alJl 199 ll) : within a radius of ~ 0.7 Mpc, its luminosity 
in the [0.5-8.0] keV bantflis ^ 1.4 x 10 45 erg s" 1 . This cluster 
exhibits several indications of a well relaxed dynamical state, for 
instance the absence of evident substructures and a central X-ray 
surface brightness peak, associated with a cool core. These charac- 
teristics support the hypothesis of hydrostatic equilibrium underly- 
ing our X-ray mass estimate. 

We performed our X-ray analysis on two datasets retrieved 
from the Chandra archive (see Table Q] for observation log); both 
the observations were telemetered in Very Faint mode, and were 
reduced accordingly. Another Chandra dataset is currently avail- 
able (ID 4974, start date: 2003-11-13, total exposure time: 58.0 ks, 

1 The X-ray luminosity was derived fitting the cluster integrated spectra in 
the radial range [0 — 0.7] Mpc using the XSPEC analysis package. 




Figure 1. True colour image of MS2137; red, green and blue colours cor- 
respond to [0.5-2.0], [2.0-4.0], [4.0-8.0] keV bands, respectively. The size 
of the field of view is 200" X 200". 



Table 1. Chandra observations summary. 



Obs.ID 


Start Time 


Total 


Net 


PI 






Expos, [ks] 


Expos, [ks] 


name 


928 


1999-11-18 08:37:37 


44.17 


20.7 


Wise 


5250 


2003-11-18 20:06:43 


41.09 


25.0 


Allen 



PI: Allen). The background light curve of this observation presents 
pronounced flares and a variable quiescent rate. A careful back- 
ground light curve screening could solve the flaring issue, allowing 
to disentangle the ICM emission from the spurious signal. How- 
ever, unlike the observations we used for our analysis, Obs.4974 
presents a widely-scattered and variable quiescent count rate. This 
feature makes it hard to distinguish the cluster emission from pos- 
sible spurious events, and could affect the background estimate. 
Moreover, the final net exposure time for this dataset was only 
~ 16 ks. So, to prevent systematic errors associated with an in- 
correct background estimate, we decided to discard this dataset. 

The data were reduced using the C1AO data analysis package 
(version 3.4) and the calibration database CALDB 3.4.2; here we 
summarize briefly the reduction procedural 
We reprocessed the level- 1 event files to include the appropriate 
gain maps and calibration products and to reduce the ACIS quies- 
cent background We used the acis.process.events tool to 
check for the presence of cosmic-ray background events, correct 
for spatial gain variations due to charge transfer inefficiency and 
re-compute the event grades. We selected the events flagged with 



2 See the CIAO analysis guides for the data reduction: 
cxc.harvard.edu/ciao/guides/. 

3 For a complete discussion on this topic, see 
cxc.harvard.edu/cal / Acis / CaLprods / vfbkgrnd / index.html. 



grades 0, 2, 3, 4 and filtered for the Good Time Intervals associated 
with the observations. We identified bright point sources using the 
script vtpdetect; the results were subsequently checked through 
visual inspection. We masked out the point sources and corrected 
for the discarded area. 

A careful screening of the b ackground light curve is necessary 
to identify background flares (Marke vitch et alj|2003l) . The back- 
ground light curve was extracted with a time bin size of ~ 1 ks in 
the energy range [2.5-7] keV, which is the most sensitive to com- 
mon flares for the S3 BI chip, where the cluster was centered in 
both observations. We applied the script analyze_ltcrv to in- 
clude only the time periods inside the 3a range. We compared the 
S3 background light curve with the one extracted in the S 1 chips 
using the energy range [2.5-6] keV, to check for undectected faint 
flares, that were not identified. The S3 background light curve was 
examined using the ChlPS facilities to identify and exclude further 
flaring events; finally we selected only time bins within a factor 1.2 
of the apparent quiescent rate, through the script lc.clean. At the 
end of the light curve screening, almost half of the initial exposure 
time was discarded (see Tab. [I). 



2.2 X-ray analysis 

Cluster mass measurements derived under the hydrostatic equilib- 
rium hypothesis are strongly dependen t on the temperature profile, 
as demonstrated bv lRasia et al.l {2006); therefore, any temperature 
bias could affect the mass estimates. We extracted the background 
spectrum from the blank-field background data sets provided by 
the ACIS calibration team, to derive ARF and RESPONSE ma- 
trices consistent both with the source and the background spec- 
trum. The blank-sky observations underwent a reduction procedure 
comparable to the one applied to the cluster data, after being re- 
projected onto the sky according to the observation aspect infor- 
mation. We derived the background re-normalization factor in the 
[9.5-1 1.5] keV band for the Observation 5250 and in the [8.0-9.8] 
keV band for the Observation 928. The normalization factor was 
derived in different bands for the two observations because of a 
difference in the energy upper limit of the FEF (FITS Embedded 
Function) calibration files associated to the observations (for the 
latter observation, the upper limit is ~ 9.8 kev). Currently, the en- 
ergy range defined for the spectrum extraction must match the one 
in the FEF files to derive a RESPONSE matrix using the CIAO 
package tool However, the observed flux in both these bands can 
be safely attributed to the background component, so the difference 
in the blank-sky renormalization bands has a negligible impact in 
the subsequent analysis. 

The background in the soft band is variable both in time and in 
space, so we verified whether the soft X-ray background derived by 
the blank-sky datasets is consistent with the observed one. If they 
were significantly different, one should take into account this factor 
in the X-ray spectral fitting. For both datasets we extracted a spec- 
trum in a region free of cluster emission, to which we subtracted a 
spectrum derived in the same region of the blank-field dataset. We 
fitted the residuals in the [0.4—1] keV energy band with a MEKAL 
model, without an absorption component and broadening the nor- 
malization fitting range to negative values. The results indicate that 
a corrective component is not necessary. 



4 See cxc.harvard.edu/ciao/bugs/mkwarf.html. 
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2.3 Spectral fitting 



We derived the cluster X-ray flux image through the combined, 
point-source removed images in the energy range [0.3-8.0] keV, 
properly dividing by the exposure maps for the two observations. 
The spectra were extracted on concentric, circular annuli centred 
on the X-ray surface brightness centroid, sized to contain at least 
5000 net source counts, up to a net count rate of ~ 40%. We se- 
lected 7 annuli, up to 2.45 arcmin. The X-ray emission appears 
well centred on the brightest cluster galaxy (BCG). The shift be- 
tween the estimated X-ray centroid and the BCG center is ~ 1 
arcsec: the uncertainty on the X-ray centroid estimate is compa- 
rable to the applied smoothing scale (1.5 arcsec). We used the 
CIAO specextract tool both to extract the source and back- 
ground spectra and to construct ancillary-response and response 
matrices. The spectra were fitted in the [0.6-7.0] keV range, ex- 
cept for the last annulus, where, due to the higher background 
level, we restricted the analysis to the range [0.6 -5.0] keV. We 
used the XSPEC software package <Arnaudlll996l) to analyse the 
spectra, that were fitted adoptin g an optically-thin plasma emis- 
sion m odel (the MEKAL model; lMewe et alJl985l ; lKaastra & Mewel 
1993), with an absorption compo nent, described b y the Tiibingen- 
Boulder model (tbabs model; IWilms et alj2"00(ih , to take into ac- 
count the effect of Galactic absorption on the observed flux. 
The Galactic abs orption was fixed to the value inferred from ra- 
dio HI maps in lDickev & Lockmanl |l990), i.e. 3.55 x 10 20 cm -2 , 
with the cluster-averaged values of Galactic absorption derived 
from the X-ray spectra fit being consistent with the radio value 
(3.9±0.8 x 10 20 crrr 2 and 4.2 ±0.9 x 10 20 cm- 2 for the Observa- 
tions 928 and 5250, respectively). Thus the free parameters in our 
model were the temperature, the metallicity and the normalization 
of the thermal spectrum. 

The spectral data were not rebinned to keep a high spectral reso- 
lution: we only grouped the channels to have at least one count in 
each channel. As a consequence, we applied the c-stat statistic to 
the spectral fitting. 

The two observations were at first analysed individually, to assess 
the consistency of the datasets and to exclude any systematic effects 
that could influence the combined analysis. A significant difference 
was detected in the spectral fit results for the most distant radial 
bin, supposedly for the higher background in Obs.5250. Thus the 
spectrum relative to Obs.5250 was discarded. We proceeded with 
the joint analysis of the two datasets, tying the temperature and 
metallicity parameters and leaving the normalization unlinked. The 
fit results from the spectral analysis of the two combined datasets 
are listed in Tab. [2] The resulting metal abundance profile presents 
a central peak with a smooth external decline, as is generally ob- 
served for dynamically relaxed clusters. Moreover, the temperature 
decrease in the central region associated with the cluster cool core 
is evident in Fig. [2] Several earlier studies stated that the presence 
of a cool core, if neglected, could be a possible source of error in 
the X-ray mass estimate, resultin g in an underestimate of the clus- 
ter virial tempera ture and mass lAllerJI 19981 ; lHallman et al]|2006l ; 
iRasia et al]|2006h . On the contrary, if the gas temperature decre- 
ment is properly taken into account, cool core clusters should be 
ideal targets for mass analyses, being typically dynamically relaxed 
systems. Since our X-ray mass estimate relies directly on the mea- 
sured gas temperature profile, the temperature decrement in the 
center is not an issue in the mass derivation. 
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Table 2. X-ray spectral fit results. The analysis was performed in the [0.6- 
7.0] keV energy range except for the last annulus, where we considered 
the range [0.6-5.0] keV. The columns list the outer radius of the annulus 
where the spectral extraction was performed, the best-fit values obtained for 
the free parameters of the model (i.e. the intra-cluster gas temperature, the 
normalization of the X-ray spectrum and the metallicity) and the c-statistic 
result (the degrees of freedom for each annulus are in square brackets). 



Spectral fit results 


Rout 


Tgas 


Normalization 


Metal abundance 


Goodness of fit 


[arcmin] 


[keV] 


xlO 4 


[2q] 




0.065 

0.114 

0.17 

0.26 

0.45 

1.0 

2.45° 


4 ,3+0.20 

4 22+ ' 19 
-0.18 

4 66+ ' 23 
*- DD -0.21 
r qc+0.28 
0.30_ Q 2 g 

5.88ip 
5.03t°;i 

„ fi o+0.59 
o - oo -0.43 


7 

' - UO -0.27 
7 £9+0.28 
' - D -0.29 
7 43+O.25 
1 ■^°-0.26 

8.67t -j 5 7 

11 17+ - 28 
11 - 1 '-0.30 

12 15+ - 32 
1Z.10_ 33 

6 n0+°- 54 


0.60±°;i? 

0.58«;» 

0.52«; 12 

50+ 12 

41+ 010 
u '^ 1 -0.09 
n oo+O.lO 

40+ 30 
u '^ u -0.21 


459.2[567] 
514.8[552] 
578. 7 [579] 
597.0[615] 
670.3[679] 
658.3[740] 
348.0[291] 



a The results obtained for the last annulus refer to a fit performed for Obs. 
928 only. 
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Figure 2. (Top panel) The X-ray spectrum of MS2137, obtained combining 
Obs. 928 and Obs. 5250. The spectrum was reduced from an aperture of 
2.45 arcmin, and shows the energy range [0.6 — 5.01 keV. We fix ed the 
Galactic absorption to the value reported in Dickey & Lockman 1 1990). The 
normalization of the model was left as a free parameter for both datasets. 
The plotted spectrum was derived with a grouping of 20 counts per bin. The 
best fit for a MEKAL model with untied normalizations is shown as a solid 
line. (Bottom panel) Fit residuals. 



2.4 X-ray mass profile 

In this work we recover the deprojected profiles of gas tempera- 
ture and density without any parametrization of the ICM proper- 
ties, which could introduce considerable sy stematic errors in the 
mass measurement results (Ra sia et al. 2006). Our estimate relies 
on the assumptions of spherical symmetry and hydrostatic equilib- 
rium and on the choice of an analytic model for the total mass. 
In the case of MS2137 these assumptions seem to be reasonable, 
since, e.g., the strongly peaked central surface brightness and a 
roughly circular X-ray emission can generally suggest a relaxed 
dynamical state. 

Here we summarize the most relev ant aspects of the technique 
we have applied; for full details, see lEttori. De Grandi & Molendil 



iMorandi, Ettori & Moscardinil d2007l) : iMorandi & Ettoril 

From the spectral analysis we derived for each annulus the pro- 
jected gas temperature T r i ng , metal abundance Z T i ng and the nor- 
malization K of the model; the latter provides an estimate for the 
Emission Integral EI — J n c n p dV ~ 0.82 J n 2 dV. To derive 
the value of these quantities in volume shells, one has to calcu- 
late the intersection of each volu me shell with 2-dimensional an- 
nuli. Geometrical considerations dKriss. Cioffi & Canizareslll983l; 
iBuotdbOOOh allow to evaluate the intersections through an upper 
triangular matrix (V) where the last pivot represents the outermost 
annulus, under the assumption of spherical symmetry. In this way 
we can recover the deprojected values applying some array opera- 
tions. 

To calculate the deprojected density profile, we combined the spec- 
tral information and the cluster X-ray surface brightness (hereafter 
SB) profile. The SB provides an estimate of the volume-counts 
emissivity _F D i n cc n 2 bin T(bin) 1//2 /DL. By comparing this ob- 
served profile to the values predicted by a thermal plasma model 
(with temperature and metallicity equal to the measured quantities 
and taking into account the effect of Galactic absorption) one can 
solve for the electron density. In this way, we obtain a gas density 
profile which is better resolved than a spectral-only profile (see the 
top panel of Fig. [3}- We computed the SB in 35 circular annuli. 
They were derived from the X-ray flux image in the energy range 
[0.5-2.0] keV, thus similarly to the spectral extraction annuli (see 
§[23}, but applying a binning of 1000 counts. The X-ray centroid 
is the same for both profiles. We rederived the mass estimate calcu- 
lating the SB profile over the range [0.3 — 8.0] keV, to verify if the 
energy range choice would affect our results. 
The total mass was then constrained as follows. The observed de- 
projected temperature profile in volume shells T s h c ii was obtained 
through the relations: 



II c 

e 

tshcii 



[(Vol T )- 1 #( J B//0.82) 

(VoFr^Lring, 

(Vol T )- 1 #(L ri n g T ring ) 

(vor)- i #(L ring ^in g ) 



1/2 



(1) 

(2) 
(3) 
(4) 



where the symbol # indicates a matrix product. 

We compare it to the predicted temperature profile T pre d, obtained 

by inverting the equation of hydrostatic equilibrium between the 

dark matter potential and the potential of the intracluster plasma, 

i.e.: 

n c M to t, mo dei(< r) _ d (n c x fcT prcd ) _ 



Gfim p - 



(5) 



r 2 dr 
where n c is the deprojected electron density. Our best-fit mass 
model was determined by comparing the observed T bs profile to 
the predicted one, which depends on the mass model parameters. 
So we minimized: 

2 



X 



E 



T. 



obs 



pred 



<y obs 



(6) 



The Tp rc d profile was rebinned to the spectral scale. In the bottom 
panel of Fig.[3]we show the projected temperature profile as deter- 
mined by the spectral analysis, the deprojected best-fit profile and 
the deprojected, rebinned profile. 

We assumed an NFW density profile, which can be expressed 

as: 

PnfwM _ &c 



(r/r s )(l+r/r s ) 2 



(7) 
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Figure 3. (Top panel) Deprojected gas density profile. Red circles mark the 
values determined by the normalization of the thermal spectra, diamonds 
represent the gas density values obtained combining both spectral and sur- 
face brightness measurements; the dashed line shows the best-fit profile. 
(Bottom panel) ICM temperature profile. Here we show the observed tem- 
perature profile (red circles, T b s in the text), directly inferred from the 
spectral analysis, the deprojected temperature values (diamonds, T pret i in 
the text), rebinned to the same intervals of the observed profile, and the 
deprojected best-fit profile (dashed line). The small scale fluctuations of 
the temperature best-fit profile reflect the local variations in the gas den- 
sity, since the condition of hydrostatic equilibrium is imposed between the 
smooth total mass profile and the gas pressure profile P gas oc n gaa X T gas 
(see Eq.[5]in the text). 



where r s is the scale radius of the halo, S c is its characteristic 
overdensity and p c is the critical density of the universe. We will 
define the NFW profile using the scale radius and the concentration 
parameter C200, i.e. the ratio of the radius where the halo matter 
density is 200 times the critical density (r2oo) to the scale radius: 
C200 = r2oo/r s . 

The best-fit values are determined by minimizing the statistical 
estimator defined in Eq. [6] A first minimization is performed in 
the ranges 0.5 < C200 < 20.0 and 10 kpc < r s < R spa t (where 
R spat ~ 676 kpc is the outer radius of the surface brightness 
profile). In this way a first estimate of the best-fit values is obtained. 
A second minimization is then performed in the Xi ± 3a range, 
being Xi the mean values inferred through the first minimization. 
The mass model parameters obtained through our X-ray anal- 
ysis, deriving the SB profile in the energy range [0.5 — 2.0] keV, are: 



Best fit X-ray results 


r s 
[kpc] 


C200 M200 R200 
[10 14 M Q ] [Mpc] 


X?od[d-o.f.] 


162.72±igj 


8.68± ;| 4.4 ±0.3 1.41 ±0.03 


2.1 [5] 



The confidence levels for the NFW parameters are shown in Fig. [5] 
The choice of the energy range over which the surface brightness 
is estimated does not significantly affect our results: for example, 
deriving the SB profile in the wider range [0.3 — 8.0] keV we 
obtain the following values: 

j+17.2 



(i) r a 

(ii) c 

(iii) M200 



153.8^20:7 kpc, 



8.93l °, 



4.1 ±0.3 x 10 14 M 



0> 



(iv) Xred[d-0.f.] = 3.1[5]. 



The gas mass profile and the total mass profile are shown in 
FigO To indicate how much the X-ray mass profile depends 
on the choice of the density profile, in Fig[4] we show also 
the mass profiles derived assuming two other parametric mod- 
els, i.e. the Rasia-Tormen-M oscardini (hereafter RTM) profile 
1 Rasia. Tormen & Moscardirnl I2004T) and the King profile ( I Kind 
19621) . The RTM model was deduced from the direct analysis of 
non-radiative A-body hydro-simulations, without any assumption 
on the gas equation of state or on the cluster dynamical equilibrium. 
The resulting DM density distribution is: 



Pn(r) = 



P0,R 



(8) 



r(r ± r*) 1 ' 5 

The King profile jKindl 19661 : fRood et al.ll 19721) was derived under 
the assumption of an isothermal particle distribution, and can be 
expressed as: 



pK(r) = 



P0,K 



(9) 



(l±(r/r c ) 2 )3/ 2 ' 

The mass estimates which were obtained adopting the RTM 
and the King mass model are M200, rtm = 5.1±0.3x 1O 14 M0 and 
M2oo,Kin g = 2.2 ±0.2 x 1O 14 M0. The corresponding x? c d values 
are 3.7 and 4.1 (over 5 degrees of freedom), respectively. When as- 
suming the RTM or the King density profile, the surface brightness 
was estimated in the energy range [0.3 — 8.0] keV. The best-fitting 
model is the NFW profile, even though the RTM profile is not a sig- 
nificantly worse fit. Clearly, the extrapolated mass derived through 
a profile shallower in the external regions is lower than the mass 
estimates obtained with more cuspy parametrizations. 

We have compared our X-ray estimate with the 
value predicted by the M — Tx relation. For instance, 
lArnaud, Pointecouteau & Prattl (2005), through the analysis 
of 10 nearby, relaxed galaxy clusters observed with XMM-Newton, 
deduced the following relations for clusters with kT 3.5 keV: 



1.49 



h(z)M 200 = 5.74 x 
h(z)M 500 = 4.10 x 



kT 
5keV ' 
kT 149 



5keV 



where h(z) is the Hubble constant normalised to its local value 
and kT is the spectroscopic temperature over the [0.1 — 0.5] R200 
region (corresponding to [0.15 — 0.75] Rsoofl The authors de- 

5 Accurate predictions on the X-ray s caling relations are also presented 
in iNagai . Kravts ov & Vikhlininl J2007l> . through the results obtained by 
IVildiluiin et alj 12003) analysing Chandra observations of 13 low redshift, 
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at z ~ 0.15. We estimate d the value of kT in a similar man- 
ner to lArnaud et al. (2005): we extracted the X-ray spectrum in 
the energy interval [0.6 — 5] keV and in the same radial range 
[0.1 - 0.5]R 2 oo. We derived the value of R 20 o = 1-37 ± 0.04 
Mpc fitting the data with an NFW profile for a density contrast 
S = 200. We considered only the Obs. 928, due to the higher back- 
ground in the external regions for Obs. 5250 (see tj |2.3t . The spec- 
troscopic gas temperature over the spatial range [0.1 — 0.5] R200 
is Tx = 4.1 ± 0.3 keV (the c— statistic value for the spectral 
fit is 293.2 over 296 degrees of freedom). This value leads to 
M 20 o,scai ^ 3.7 x 10 14 M Q and M 50 o,scai 2.6 x 10 14 M Q . These 
estimates are lower than ours (i.e. M200 ^ 4.4 ± 0.3 x 1O 14 M 
and M500 ^ 3.5 ± 0.2 x 10 14 M Q ), but are compatible at a 4cr level 
with our results. The X-ray mass of MS2137 seems to be consistent 
with the M—Tx scaling relation found for relaxed galaxy clusters. 



Figure 4. MS2137 total mass profiles, as determined through the X-ray 
analysis for different choices of parametrised mass models. The vertical 
line indicates the boundary of the spectral extraction, the filled circles mark 
the value of R200 for the corresponding mass profile. The measured gas 
mass profile is also shown as a green solid line. 



O 
Oh 
M 



«3 
«3 




200 



Figure 5. Confidence levels determined through the X-ray analysis for the 
NFW profile parameters C200 and r s . The contours indicate the la, 2cr and 
3<r c.f, and the black circle marks the best-fit values. 



rived the temperature Tx from a single-temperature fit to the in- 
tegrated spectrum. The highest redshift cluster in the sample is 



relaxed clusters. The scaling relations in iNagai et alj 120071) are calibrated 
with the gas temperature derived in the wider r adial range 70 kpc < r < 
R500 . We applied the scaling relations found bv lArnaud et al l j2005l) to re- 
fer to a spatial ran ge where our s pectro scopic temperature estimate is more 
robust. However, lArnaud et alj 120051) asserte d that their s caling relation 
estimates are very similar to those derived in INagai et alj j2007l) . so our 
conclusions would be unchanged. 



3 STRONG LENSING ANALYSIS 
3.1 Overview 

MS2137. 3-2353 is a well-known lensing system; the rad ial feature 
in its double-arc configuration is the first ever detected dFort et all 
1992). The radial arc may be very useful in constraining the cluster 
central density slope. For axially simmetric lenses, the radial 
and tangential critical curves arise if A r = 1 — (d/dx)(m/x) 
and A t = 1 — m/x 2 vanish, respectively, A r and A t being the 
eigenvalues of the Jacobian matrix of the lens mapping (see 
Bartelman n & Naravanlll995l for an extended review). Thus, the 
tangential and radial arc positions allow one to constrain the en- 
closed total mass and the mass distribution derivative, respectively. 
We carried out a strong lensing analysis on deep HST data 
retrieved from the Space Telescope Archival The data consists 
of a WFPC2 association image, taken with the F702W filter and 
obtained combining 10 exposures for a total integration time of 
22.2 kfl (a portion of the HST image is shown in Fig.[6j. 
The tangential and radial arc distances from the BCG center are 
~ 15 and ~ 5 arcsec: the former is composed of two images with 
reverse parity ( API and A02 in Fig. [6] following the notation of 
iGavazzil l2005bl) . associated with the counter-images A2 and A3. 
A fifth demagnified image lying near the center of the BCG is 
predicted if the central density slope is shallower than the isother- 
mal model slope. Thus the position of the central image would 
provide an additional c onstraint on the inner density (for example, 
see the results of iGavazzi. Fort. Mellier. Pello & Dantel-Forj 
2003), and as a consequence on the BCG mass. However, the 
BCG-subtracted image presents residuals in the central region, 
where the fifth image, if present, is predicted. The noise level after 
the BCG subtraction makes a clear identification of the central 
image difficult. For these reasons the central counter-image was 

not included as an additional co nstraint. 

Spectroscopic measurements dSand, Treu & Ellisl l2002h stated 
that the two groups of images arise from 2 different sources 
at redshift gtang = 1.501 a nd 2 ra( j = 1.502; recent data 
dSoucail. Covone & Kneibll2007l) confirmed that the counter-image 
of the radial arc Bl is the feature B2 (see Fig.[6]l. 
In the next section we will present a parametric strong lensing 
reconstruction of the cluster mass distribution, performed using the 



6 archive.stsci.edu/ 

7 Proposal ID: 5402; PI: I.Gioia. 
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Figure 6. Multiple image interpretation of the lens system, overplotted on 
the HSTAVFPC2 image of MS2137 (filter: F702W; number of combined 
images: 10; total integration time: 22.2 ks; proposal ID: 5402; PI: I.Gioia). 
The size of the field of view is ~ 42 arcsec X 45 arcsec; the coordinates 
of the images used as constraints in the strong lensing analysis are listed in 
Tabled 



Table 3. Coordinates of the multiple images considered in the MS2137 
strong lensing modeling. Coordinates are in arcsec with respect to the BCG 
center. We did not include as constraint the fifth central imag e, that was not 
clearl y identified. Spectroscopic measurements jSand, Treu. Smith & Eltisl 
120041) indicate a redshift of ztang — 1.501 and z rac i = 1.502 for the first 
and the second system, respectively. 

Multiple-image systems 



Id 


Im.l 


Im.2 


Im.3 


Im.4 


l.a 


-5.20 


13.87 


2.61 


14.86 


13.59 


-1.08 


-12.03 


-15.29 


l.b 


-4.88 


13.98 


2.31 


14.89 


13.46 


-1.26 


-12.36 


-14.94 


l.c 


-7.32 


13.68 


0.04 


15.77 


13.06 


-2.96 


-13.32 


-13.63 


l.d 


-7.45 


12.74 


4.69 


14.80 


13.62 


-2.13 


-11.83 


-15.00 


l.e 


-3.79 


14.52 


0.27 


15.06 


13.33 


-1.40 


-12.62 


-14.71 


l.f 


-6.85 


13.86 


-1.35 


15.71 


12.93 


-2.78 


-13.51 


-13.68 


1-g 


-7.81 


13.26 


1.32 


15.82 


12.86 


-3.28 


-13.32 


-13.50 


l.h 


-5.70 


14.55 


-2.14 


15.60 


12.80 


-3.06 


-13.55 


-13.53 


Li 










13.13 


-1.23 


-13.28 


-14.45 


lj 










12.01 


-2.12 


-15.13 


-12.80 


l.k 










-14.17 


-13.40 


12.50 


-2.01 


2.a 


0.04 


6.79 


0.33 


3.41 


-7.28 


-22.49 






2.b 


0.00 


6.44 


0.27 


3.88 


-7.22 


-22.75 






2.c 


-0.06 


5.87 


0.17 


4.21 


-7.05 


-22.88 






2.d 


-0.03 


5.28 


0.10 


4.50 


-6.92 


-23.04 







Lenstoo jf| analysis software dKneibet all 1996; Juli o et al]2007l) . 



3.2 Lens modeling 

To identify the multiple ima ge systems, we st arted by referring to 
the configuration provided bv lGavazzll J2005bh (hereafter G05), the 



The code is available at the url www.oamp.fr/cosmology/lenstool. 



most detailed knots catalogue available in the literature. We then 
iteratively refined the system by selecting several sets of ellipses, 
likely belonging to a common source area (their coordinates are 
listed in TableO. 

Our lens mass model consists of a cluster-scale mass component 
and of two galaxy-scale components, namely the BCG and the 
galaxy near the tangential arc. The latter (flagged 01 in Fig. [6]l 
was included to account for its li kely perturbative effec t, due to its 
proximity to the tangential arc dMeneghetti et al. 2007). We mod- 
eled the cluster scale component using an elliptical NFW profile, 
as implementend in the Lenstool code (for a det ailed discussion on 
the lens properties of this density profile, see Ib artelmann: 1 1 996J : 
Meneghetti, B a rtelmann & Moscardinil fetXBbl : for full details on 
the e lliptical NFW implem entation in Lenstool the reader can re- 
fer to lGolse & Kneibll2002h . 

Our main clump model has seven free parameters: the coordinates 
of the center of mass (x c , y c ), the ellipticity of the mass distribution 
(e), the position angle (6), the characteristic velocity dispersion oq 
and the NFW profile parameters (the concentration C200 and the 
scale radius r s ). The center of mass was initially set to the BCG 
centroid, but it was allowed to vary by ±15 arcsec; the ellipticity 
and the position angle of the dark matter clump were also, as a first 
estimate, set equal to the BCG values, then were optimized in the 
ranges [0.0 < e < 0.5] and [130.0 sC 9 sC 160.0] degrees. The op- 
timization ranges for the scale radius and the concentration param- 
eter are [50 ^ r s ^ 650] kpc and [0.5 ^ C200 ^ 20.0]. The choice 
of these intervals is motivated by cluster analysis results from pre- 
vious studies. Of course, strong lensing allows to constrain only 
the inner matter distribution, so the estimate of the virial radius - 
and therefore of the concentration parameter - using strong lensing 
alone is less robust; however, we decided not to fix any mass distri- 
bution parameter, but rather to sample them in wide ranges to avoid 
biasing our results. 

The galaxies were modeled either as dual Pseudo Isothermal El- 
liptical Mass Distribution - hereafter dPIE - with a null core ra- 
dius, _oj^_as_erfij3tical^ngjfJ^^^ hereafter SIE 



(see Nataraian & Kneibfl997l . lLimousin. Kneib & NataraiarfeOOSi 
Elfasd ottir et alj|2007l for the properties of these potentials and for 
full details on th eir implementati o n in the Lenstool code). Un- 
like the studies of lGavazzil J2005bl) ; [Sand et alj J2008), we did not 
model the BCG stellar velocity distribution in our analysis. We as- 
sumed a parametric model for the BCG mass (dPIE or SIE), so 
its mass is determined by the central velocity dispersion (and the 
scale/core radius values in the dPIE case). We are aware that stellar 
dynamics data would help to solve the existing degeneracy between 
the stellar mass and the mass of DM and gas in the cluster center. 
However, here we are more interested in determining the shape of 
the total mass profile at radii larger than those where the stellar 
mass is a non negligible mass component. 

We fixed the galaxies ellipticity and orientation to the 
values derived by the HST image, using the IRAF external 
package isophote of the stsdas distribution. When mod- 
eling the galaxies as dPIE mass distribution, the values we 
imposed to the dPIE scale radius are almost the same of 
ISand. Treu, Ellis. Smith & Kneibl ( 120081) . so r cut ,BCG= 22 kpc and 
r C ut,gal= 5 kpc. The other free parameter in the galaxy optimiza- 
tion was the central ve locity dispersion op ; the sp ectroscopic mea- 
surements obtained by San d. Treu & Ellis] |2002) with the Echelle 
Spectrograph and Imager on the Keck II telescope indicate a cen- 
tral value for the BCG stellar velocity of cr mea sured,BCG — 325 ± 
35 km s _1 . However, the cD measured velocity dispersion cannot 
be directly related to the "characteristic" one in galaxy mass mod- 
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Figure 7. Images predicted by our best-fit lens model. In this case the 
galaxies included in the fit (i.e. the BCG and the galaxy flagged 01 in 
Fig. |6) were modeled as dPIE potentials. The size of the image is ~ 
42 arcsec X 45 arcsec. The red [yellow] lines represent the critical [caustic] 
lines. 

els, since the observed dispersion depends upon the unknown dark 
matter potential contribution to the stellar dynamics. The uncer- 
tainty on the velocity dispersion parameter could poten tially intro - 
duce additional errors in the strong lensing modelling (Wu 1993). 
Despi te of this, the BCG stellar mass is an importan t lens compo- 
nent (Meneghett i. Bartelmann & Moscardinil l2003al iHilbert et al.l 
2008), especially when modeling radial arcs, which, due to their 
position very close to the cluster center, are extremely sensitive 
to the cuspy distribution of lu minous matter (see the earlier find- 
ings of Miralda-Escude 1995 on this topic). We thus imposed the 
measured cr mcasurc d,BCG as an upper limit for the velocity dis- 
persion optimization range, and set as lower limit the value of 
0"o, bcg — 230 km s -1 , that is plausible given the high luminosity 
of the BCG. A discussion on the effect of this choice on the lensing 
results is presented in § 13.31 

The optimization procedure of the Lenstool code is performed in 
the source plane, by mapping the positions of the input multiple 
images back to the sour ce plane, and requiring they have a minimal 
scatter jjullo et al.ll2007l) . The goodness-of-fit is evaluated through 
the following \ 2 estimator: 

3 

where the subscript j refers to the multiply imaged systems. 

We assumed the uncertainty in the lensed image positions to 
be a i — 0.3": mapping it back to the source plane one obtains an 
estimate of the error in the source plane. 

3.3 Strong lensing results 

The optimization and the error measurement are performed through 
the Bayesian optimization method implemented in Lenstool, im- 
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Figure 8. Confidence levels for the NFW mass model parameters deter- 
mined through a strong lensing analysis. The contours show the parameter 
non-marginalized probability density, and refer to the model [1] in Table 
\4\ [Left panel] NFW concentration parameter C200 vs the scale radius r s ; 
[right panel] the concentration parameter vs the cluster total mass M200, 
extrapolated up to R200 ■ 



Table 5. Best fit parameters for the galaxies included in the lens system. The 
galaxies were parameterized as dPIE or SIE potentials in the lens model [1] 
and [2], respectively (in Tab.|4]we list the fit results for the other parameters 
of the lens model). We report the mode of the distributions, with the 68% 
statistical errors, and the best-fit values enclosed in round brackets. Values 
in square brackets were not optimized. 



Galaxy model parameters 



Model 


Halo 


Ellipticity 


Rcut 

[kpc] 


CO 

[km s _1 ] 


[1] 

[2] 


BCG 
Galaxy 

BCG 
Galaxy 


[0.18] 
[0.38] 
[0.18] 
[0.38] 


[22.0] 
[5.0] 


230.7^1(230.5) 
150.2t23;| (159.2) 
231.4± 33 3 (240.3) 
124.5±}g;5(124.2) 



posing an optimization rate of S\ = 0.1. 

A summary of the best-fit values for two lens configurations is re- 
ported in Table [4] 

In Table [5] we list the galaxy model values inferred through 
the strong lensing optimization. The configuration we will take as 
a reference in the following is labelled [1]. 

The lensed images predicted by the reference model are shown in 
Fig(7]and the corresponding likelihood contour levels for the NFW 
parameters are shown in Fig[8] In Table [6] are reported the \ 2 an d 
the mean scatter for each image system. 

The curvature of the radial arc, peculiar of this strong lensing 
system, is not well reproduced by our lens models. Due to its po- 
sition very close to the cluster center, the radial arc is affected by 
the BCG mass and dynamics. Thus a possible explanation of the 
offset between the predicted and the observed radial arc curvature 
could be, e.g., the effect of a potential with non-constant ellipticity 
or position angle, that cannot be take n into account with the cur- 
rent mass model implementations (see lGavazzi et alj |2003). Other 
features of the arc (like its position or its radial extension) are quite 
well reproduced by our models. So, despite of the uncertainty on 
its curved shape, the radial arc provides indeed a tight constraint on 
the position of the radial critical line, and as a consequence on the 
total mass derivative in the cluster central region. 

Resuming the strong lensing results, the best-fit models indi- 
cate a negligible shift of the cluster center with respect to the BCG 
centroid, and a slight offset between the BCG and the cluster halo 
projected orientation , A8 = 5°, which is lower than the AO = 13° 
found by G05 and bv lSand et al] |2008). Such a misalignment could 
have either a geometrical explanation (e.g. the projection effect 
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Table 4. Best fit parameters, obtained for different lens models. Hereafter, we will refer to the first model as the reference model. Quoted errors correspond to 
lcr confidence level. The columns indicate the lens configuration (the cluster halo model + the galaxy model), the coordinates of the center of mass, expressed in 
arcsec with respect to the BCG centroid, the ellipticity and position angle (WCS aligned) of the cluster clump, the scale radius and the concentration parameter 
(both referred to a halo overdensity of 200), the extrapolated mass up to R200 > the characteristic velocity dispersion of the NFW profile implementation in 
Lenstool and the total x 2 . calculated over 15 knots systems, for a total of 50 image areas. When the best-fit value is significantly different from the mode of 
the distribution, the former is reported enclosed in brackets. The errors on x c , y c and on the ellipticity are ~ 0.1 and ~ 0.01, respectively. Ellipticity e here is 
expressed as (a 2 — b 2 )/(a 2 +b 2 ). 



Best fit strong lensing results 



Model 


X c 

[arcsec] 


[arcsec] 


e e 

[deg] 


r s 
[kpc] 




C200 M200 CT 

[10 14 M©] [kms- 1 ] 


X 2 ot [d.o.f.] 


[1] NFW+dPIE 


0.0 


0.0 


0.10 147.2±g;|| [147.4] 148.9; 


(-22.0 
-24.7 


[143.5] 9.50; 


^;f 7 [9.76] 4.1lt°;2>[4.31] 1377.3±||;| [1366.4] 


17.6[61] 


[2] NFW+SIE 


0.1 


-0.1 


0.10 149.3 ± 0.7[149.2] 118.3; 


(-26.0 
-14.3 


[117.3] 10.47; 


tl;lg [10.98] 3.38l ;i[3.36] 1291.2l 28 j! 


10.9[61] 



Table 6. The results of the source plane optimization performed with the 
Lenstool code for our reference model, labelled [1] in Table[4] The columns 
list the knots system, the \ 2 of the system, the root-mean-square computed 
in the source plane (rms s ) and the corresponding one in the image plane 
(mis;). The last row reports the total x 2 estimator, referring to a model 
with 61 degrees of freedom. 



Goodness of fit 



System 


x 2 


rms s ["] 


rmsi["] 


La 


0.48 


0.034 


0.27 


L.b 


1.72 


0.064 


0.37 


l.c 


2.18 


0.075 


0.22 


l.d 


0.90 


0.053 


0.23 


l.e 


1.80 


0.062 


0.50 


L.f 


1.61 


0.062 


0.86 


l-g 


1.06 


0.057 


0.37 


L.h 


0.21 


0.021 


0.25 


l.i 


0.11 


0.030 


0.08 


1.1 


0.16 


0.038 


0.10 


l.m 


0.14 


0.034 


0.07 


2.a 


3.40 


0.168 


0.33 


2.b 


1.97 


0.114 


0.35 


2.c 


1.13 


0.078 


0.28 


2.d 


0.75 


0.057 


0.44 


Total 


X 2 W.o.f] 


rms s ["] 


rmsi["] 




17.62 [61] 


0.072 


0.39 



325 km/s. Under these conditions we derived the following values: 



Parameter 


Unit 


Mode 


Best fit 


x c 


arcsec 


0.2 ±0.1 




y c 


arcsec 


0.2 ±0.1 




e 


l's 

C200 

M200 

O"0,halo 
CO, BCG 


deg 
kpc 

1O 14 M 
km s~ x 
km s~ x 


0.11 ±0.01 
146.21°;! 
153.1+ 2 ^ 

9.2it8:S 

4.40l° ;« 
1359.0l 3 2 ^ 
327.0lg; 2 5 


[160.1] 
[8.93] 
[4.59] 

[325.4] 


X? ot [d.o.f.] 




61.1[61] 





The values in square brackets indicate the best-fit results if 
they are significantly different from the mode of the distribution. 
The parameter values in this case are within 2a respect to the ones 
presented in Tab.|4] The results obtained decreasing the BCG veloc- 
ity dispersion (and thus its mass) are shown in §[5] So, considering 
the results obtained for the two boundary conditions (no halo con- 
tribution to the BCG potential well and no mass in the BCG) we 
can estimate that, regard to the strong lensing mass, the systematic 
errors associated with the uncertainty in the BCG mass account for 
an additional ~ 25% in the final error budget. We would like to 
stress again that this additional error component arises in the case 
of galaxy clusters with central lensed features, when trying to dis- 
entangle the central galaxy mass from the underlying halo mass 
and when no additional constraint (e.g. the BCG stellar dynamics) 
is included. 

Anyway the NFW parameters are quite well constrained by the 
cluster strong lensing features, and we did not detect a high concen- 
trati on for this cluster, as stated i n previous works in the literature 
fsee lComerford & Nat araian 2007 and §|5]for further discussion). 



of a prolate halo -see G05- or the difference of halos triaxialities 
dKeeton. Kochanek & SeiiaSl 19971) ) or a dynamical origin, linked 
to the B CG formation process (|C/uadri. Moller & Nataraian 2003; 
see also lMinor & Kapl ingha 3120071) . 

To estimate the magnitude of the systematic error on the lensing 
mass introduced by the uncertainty on the velocity dispersion, we 
optimized again the lens model imposing as lower limit <to,bcg = 



4 X-RAY AND STRONG LENSING COMPARISON 

Coupling the X-ray and lensing mass estimates is not a trivial task, 
as already stated by G05, due to their different dependencies on the 
characteristics of galaxy clusters (like asphericity) as well as for 
the still incomplete treatment of the biases that could affect both 
measurements. Though joint analyses of homogeneous effects (e.g. 
weak and strong lensing) often supply converging estimates, lead- 
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Table 7. Best fit parameters for the NFW profile. The strong lensing results refer to the reference model, labeled [1] in Tab. [4] Quoted errors correspond to 
lcr confidence level. For a comparison, we reported the r esults obtained in s om e earlier representative st udies: the original values were converted according to 
common definitions. We list here the findings obtained by Allen et al. (2001) and Schmidt & Allen 1 2007) through the X-ray analysis of Chandra data (columns 
A01 and S07, respectively) and by Gavazzi l2005b|) and lComerford et alj d2006t) through a strong lensing study (columns G05 and C06, respectively). Here, 
again, e = (a? - b 2 )j(a 2 + b 2 ). 

Comparison with previous results 





X-RAY 


SL 


A01 


S07 


G05 


C06 


Units 


r s 

C200 

M 200 


162.72+^J? 

8 .68i8:S 

4.43 ± 0.25 


148.9+24 ? 

n -0+1-29 
y - ou -0.87 
4 11+0.75 


160.0 ± 30 

°- a -0.92 


180.0 ± 20.0 

o 1Q +0.54 


158.0+i| 

11 QO+ - 77 

7 5fi+ - 63 
'■°°-0.54 


91.4 ±3 
13.0 ±1.0 
2.9 ±0.4 


kpc 
1O 14 M 


e 




0.1 ±0.01 






0.25 


0.12 





ing to remarkable results feradac et aT. 2005, 200ol lLirnousin et"aL1 
2007; iMerten et alj|2008l) . comparing 2- and 3-dimensional mea- 
surements requires further caution. Recently, several interesting 
attempts have been made to infe r cluster masses thro ugh multi- 
wavelength, combined analysis ( Mah davi et ail I2007T) . in some 
case s suggesting that a direct comparison needs further investiga- 
tion dLemze et alj|20081) . We decided instead to perform two paral- 
lel analyses to understand if these methods lead to discrepant mass 
estimates forMS2137. 

We found that our X-ray and strong lensing results are in good 
agreement: both the extrapolated value of the total mass M200 and 
the mass model parameters agree within the 2a range. The param- 
eter probability distributions are mutually consistent (see the left 
panel of Figure [5J, although the concentration inferred from the 
lensing analysis is slightly higher than the X-ray one. A mild elon- 
gation of the cluster along the line of sight could be a possible ex- 
planation, as already suggested by G05. Another reason could be 
the additional uncertainty associated to the BCG mass. To evaluate 
the differences between the X-ray and the strong lensing results, we 
show in the right panel of Figure|9]the projected total mass enclosed 
in concentric cylinders inferred from both analyses. The assessment 
of the errors on the strong lensing mass was performed by stacking 
500 mass map realizations, obtained through the Lenstool Bayesian 
optimization, and computing the standard deviation at each map 
point. 



5 COMPARISON WITH PREVIOUS STUDIES 

In this section we will briefly compare our resul ts with some pre- 
vious studies, namely the findings obtained by 1 Schmidt & Allenl 
d2007l) (hereafter S07) and I Allen. Schmidt & Fabianl d200ll) (here- 
after A01) through the X-ray analysis of Chandra data, and 
those inferred t h rough a strong lensing analysis by G05 and 
Comerford et alj d2006l) (hereafter C06: see also iGavazzi et all 
20031 : IComerford & Natara^l200llshu et al.ll2008h . 
S07 studied the same two Chandra datasets considered in this 
work, whereas the analysis of A01 was carried out on the only 
Chandra dataset available at that epoch (Obs. 928). The method 
applied in their studies are similar to ours. Briefly, S07 and A01 
derived the observed X-ray surface brightness profile and the de- 
projected X-ray gas temperature profile, and combined them to 
determine the gas mass and total mass profile of the galaxy clus- 
ter, under the assumptions of hydrostatic equilibrium and spheri- 



cal Symmetrel (for more details, s ee also Iwhite. Jones & Format] 
1 19971: lAllen. Ettori & Fabianll200l|) . The gas temperature and den- 
sity profiles for the single case of MS2137 were not listed in S07 
and A01, since in both cases the authors presented the results of 
careful analyses of large samples of objects. Thus a direct compar- 
ison of the derived profiles was not possible. Concerning the final 
estimates of the NFW parameters, there is a good matching be- 
tween our results and the values presented both in S07 and A01. 
Despite some differences in the analyses (for example the spectral 
binning and the X-ray background estimate) the results are mutu- 
ally consistent within a 2a level (see Tab(7}. 
A comparison with the previous strong lensing analyses in- 
volves an extended d i scussi on, since, as already stated by 
Comerford & Nataraia nl d2007l) . the scatter in the differing results 
found in the literature regarding this cluster is quite large. A com- 
plete discussion on this topic is beyond the aim of this paper: here 
we will simply compare the results and the methods of analysis, 
and possibly suggest some explanations. For full details on the ba- 
sis of t he optimization m eth ods, the reader can refer to the related 
papers, iGavazzil d2005bh and lComerford et alj d2006l) respectively. 
Both the analyses of G05 and C06 were performed exploiting their 
own inversion codes. The source plane \ 2 estimator adop ted in 
G05 fo llows the Xs definition in the lensmodel software dKeetonl 
l2001allbh . G05 reported a value for the scale radius r s consistent 
with our estimate, whereas there is disagreement between our esti- 
mates of the halo concentration and mass. We couldn't find a con- 
vincing explanation for the difference in our results. 
The code utilised in C06 optimizes the lens model parameters by 
minimizing a figure-of-merit function, defined as the sum of three 
components both in the image and in the source planes. The op- 
timization procedure is meant to shrink: (i) the distance between 
each data point position and the nearest predicted one, (ii) the dis- 
tance between the predicted point positions and the nearest data 
point (thus penalising models that produce unobserved images) and 
(iii) the size and n oncompactness of the predicted sources (see C06 
and lShu et alj|2008l for a detailed discussion). 

This fitting technique is very accurate: however, in their anal- 
ysis C06 modeled the lens as a single elliptical NFW potential, 
not including any galaxy component. Thus, the free parameters of 
their lens model are the scale radius r s , the ellipticity, the posi- 
tion angle 6 and the scale convergence k s = p s r s /E cr it, where 



9 S07 carried out two separate mass analyses, deriving a mass estimate of 
the cluster total matter and of its three components (dark matter, gas and 
cD galaxy). In the following, we will refer to the results obtained in the first 
case. 
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Figure 9. [Left panel] 3<r marginalized probability distribution for the NFW parameters - concentration and scale radius. The regions refer to the the X-ray 
(drawn in green) and strong lensing (drawn in red) results. The square [triangle] indicates the X-ray [strong lensing] best fit values. [Right panel] Projected 
total mass enclosed in cylinders of radius R. The solid green line represents the projected total mass profile derived through the X-ray analysis (the grey area 
indicates the lcr region), the red line with error bars show the result of the strong lensing analysis. The error bars indicate the lcr errors inferred from the 
Bayesian optmization. The green [red] dashed horizontal line indicates the spatial range over which the X-ray [strong lensing] analysis has been performed. 



as usual E cr it represents the critical surface mass density. The re- 
sulting MS2137 mass parameters are r s ~ 91.4 ± 3 kpc and 
c ~ 13.0 ± 1.0, leading to a total mass extrapolated up to R200 
of M200 ^ 2.9 ± 0.4 x 10 14 M Q . Thus, the total mass and scale 
radius estimates of C06 are lower and their inferred concentration 
parameter is higher than our best-fit values. 

We perf ormed a single m ass component analysis, in analogy with 
C06 and lShuetall d2008l) . The best-fit values we obtained in this 
case are: 

(i) r a = 92.8 ± 10.0 kpc, 

(ii) c = 13.5 ± 1.1, 

(iii) M200 = 3.1 ± 0.3 x 10 14 M®, 

which are in excellent agreement with the C06 values; the to- 
tal x 2 is 29.6 for 63 d.o.f. This result highlights once again that 
the BCG has a n on-negligible effect on strong lensing properties 
(Hilbert et al. 2008) when dealing with lensed features sensitive to 
the central mass distribution, such as radial arcs. Ignoring the BCG 
stellar mass component can return an "artificial" high-concentrated 
profile, to account for the missing central galactic component. 



6 DISCUSSIONS AND CONCLUSIONS 

We have presented new X-ray and strong lensing mass measure- 
ments for the galaxy cluster MS2137. 3-2353, taking advantages 
of the high-resolution Chandra and HST data. Our X-ray analy- 
sis benefits from a non-parametric study of the gas temperature and 
density profiles. By combining them, we have recovered the cluster 
mass profile adopting the NFW model, under the assumptions of 
spherical symmetry and hydrostatic equilibrium. 
The strong lensing reconstruc tion was performed through the pub - 
licly available code Lenstool dKneib et aP l 19961 : ljullo et alj[2007h . 
which allowed for a parametric analysis comparable to our X-ray 
study. Our mass estimates are consistent within the errors, lead- 
ing to a mean value of the total mass, extrapolated up to R200, 
of ~ 4.4 ± 0.3 x 1O 14 M0. We did not find a strong discrepancy 
between the strong lensing and the X-ray NFW parameter results, 
since the best-fit values mutually agree at a 3er confidence level. 
The probability distributions are mutually compatible. The slight 



difference in the allowed parameter space could be explained, for 
example, by a mild cluster elongation along the line of sight, which 
would affect mainly the strong lensing results. Moreover, the uncer- 
tainty on the BCG mass determination could introduce a systematic 
bias in the lensing measurements that we could quantify as an addi- 
tional uncertainty of ~ 25%, since the extrapolated strong lensing 
mass estimate in this case depends upon the BCG stellar mass bud- 
get. 

The agreement between these two mass estimates for a relaxed 
cluster, like MS2137, supports the statement that, for this kind of 
system, both the strong lensing and the X-ray analyses provide an 
unbiased way to estimate the cluster mass. As a general observation 
on the physics of this cluster, the convergence of these mass results 
implies that its main pressure support mechanism against gravity 
is the thermal pressure, so the assumption of hydrostatic equilib- 
rium is verified. Non-thermal processes, such as turbulent motions 
or magnetic field effects, seemingly have little impact on the equi- 
librium state for this cluster, that appears a good example of relaxed 
cluster prototype. To verify if MS2137 represents a peculiar case in 
the cluster mass discrepancy debate, or if the relaxed dynamical 
state constitutes a sufficient condition for the convergence of the 
strong lensing and X-ray mass estimates, further investigations of 
both single cluster cases and large uniform samples are fundamen- 
tal, and would provide new, useful insights on the galaxy cluster 
physics. 
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